Ultra-High Energy Cosmic Rays and the GeV-TeV Diffuse Gamma-Ray Flux 
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Ultra-high energy cosmic ray protons accelerated in astrophysical objects produce secondary elec- 
tromagnetic cascades during propagation in the cosmic microwave and infrared backgrounds. We 
show that such cascades can contribute between ~1% and ~50% of the GeV-TeV diffuse photon flux 
measured by the EGRET experiment. The GLAST satellite should have a good chance to discover 
this flux. 



■ 

o ■ 
o : 

(N . 

Oh: 
<: 



43 ■ 

Or 
i ' 

o ■ 



> 

(N 
O 

o 



X 



Introduction. Recently the HiRes collaboration estab- 
lished [l[ the existence of the GZK cutoff 0]. If con- 
firmed, this result suggests an astrophysical origin of 
ultra- high energy cosmic rays (UHECR). 

There are two important contributions to secondary 
electromagnetic (EM) cascades from UHECR. One comes 
from the GZK process of pion production in interactions 
of UHECR protons with cosmic microwave (CMB) pho- 
tons. Most of the energy transferred to photons, electrons 
and positrons in the subsequent pion decays would cas- 
cade down to GeV-TeV energies, at which the Universe 
is transparent to photons. If the spectrum of primary 
protons is a power law oc E~ a , pion production energy 
losses and thus the energy deposited into EM cascades 
increases with decreasing power law index a. 

The second source of EM cascades is pair production 
by protons on low energy photons, p + 7 — > p + e + + eT . 
This process is more efficient for steeper injection spectra 
with larger a. Both processes together imply a minimal 
secondary photon flux at GeV-TeV energies. 

The goal of the present work is to study the con- 
tribution of EM cascades from UHECR proton inter- 
actions with background photons to the diffuse 7— ray 
background in astrophysical scenarios. We study the pa- 
rameter space of UHECR models which fit the HiRes 
energy spectrum with the GZK cutoff. We show that 
the UHECR contribution to the 7— ray flux is in the 
range 1-50 % of the diffuse flux measured by EGRET 0. 
Relatively high values of this flux should enable the 
GLAST satellite [5] to disentangle it from other contribu- 
tions such as from starforming galaxies Q .starbursts 7 1 , 
large scale structure formation shocks [8|, AGNs 9|, 
blazars [l(J, and 7— ray bursts 11 1. 

Modeling the primary proton flux. We parametrize the 
proton injection spectrum as dN/dE oc E~ a 0(E max — 
E), where E max is the maximal proton energy and a is 
the power law index for which we consider the ranges 
2 x 10 20 eV < E max < 10 21 eV and 2 < a < 2.7, 
respectively. 

We also consider the case of sources with variable den- 
sity and/or luminosity. We assume the comoving source 
density to scale as n(z) = no(l+z) m 9(z niax —z)8(z — z min ) 
where m parameterizes the luminosity evolution. We 
consider the range — 2 < m < 4, which practically in- 



cludes all astrophysical scenarios. The parameters z m i n 
and z max are the redshifts of the closest and most distant 
sources, respectively. We choose z min < 0.01 to avoid a 
GZK cutoff more pronounced than observed [l| . We fixed 
z max = 3 which is large enough to take into account cos- 
mologically distant sources. We note that both the lu- 
minosity evolution index m and the maximum redshift 
£max have similar qualitative influence on both proton 
spectrum and GeV-TeV cascade fluxes. 

For the rectilinear propagation of protons and cascades 
of secondary electrons, positrons and photons we used 
two independent codes [TH, [l3j], which we compared on 
the level of individual interactions. For pion production 
by protons and neutrons both codes use the SOPHIA 
generator [3]. In addition, protons loose energy due to 
production of e + e~ pairs, while neutrons decay. Sec- 
ondary photons produce single and double pairs on the 
low energy photon background. Electrons and positrons 
interact via inverse Compton scattering and triplet pair 
production, and undergo synchrotron energy losses in ex- 
tragalactic magnetic fields. Due to all those reactions sec- 
ondary photons, electrons and positrons cascade down to 
TeV energies, where they are not affected by interactions 
with the CMB any more. However, even at these energies 
they can still interact with infrared (IR) and optical pho- 
ton backgrounds. We use the recent model of Ref. [HI 
for these backgrounds. 

Highly structured sources and large scale magnetic 
fields can lead to enhanced synchrotron fluxes up to TeV 
energies [l6j . To be conservative we neglect magnetic 
fields and 3 dimensional effects here. 

We fit the HiRes spectrum with the method described 
in Ref. [17J. Among all models characterized by m, a, 
£ max , z m in , and Zmax we choose those which fit the latest 
HiRes spectrum [l| at the 95 % confidence level, taking 
into account empty bins above the highest energy events 
observed, as well as an energy uncertainty AE/E = 17 
%, which influences the shape of the spectrum around 
the GZK cutoff. 

We consider the two main scenarios for the transition 
from a cosmic ray flux dominated by galactic sources to 
one dominated by extragalactic sources: In the "dip sce- 
nario" extragalactic protons dominate over the galactic 
contribution down to a few 10 17 eV and the dip observed 
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in the energy spectrum between ~ 1 and ~ lOEeV is 
caused by pair production by these protons We re- 
quire predicted fluxes in this scenario fit the HiRes spec- 
trum for E > 2 EeV. 

In the second scenario, the "ankle" in the spectrum 
around ~ 5 x 10 18 eV is due to a cross-over from low- 
energy galactic to high-energy extragalactic cosmic rays. 
Recent versions of this model include a mixed composi- 
tion of UHECR at the highest energies [HI . To be con- 
servative, in this case we fit the HiRes spectrum only for 
E > 40 EeV. 
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FIG. 1: Primary proton and secondary 7—ray fluxes for an 
injection spectrum oc E~ 2Ar> up to 10 21 eV (shown as red, 
solid straight line marked Pmj) which evolves as (1 + z) 3 be- 
tween z m j n = and z max = 3. In this case the UHECR 
flux is dominated by extragalactic protons down to a few 
10 17 eV The lower red (solid) line is the proton flux, the 
blue (dash-dotted) line is the corresponding secondary 7—ray 
flux. The green (dashed) line is the proton flux without e 
production and the magenta (dotted) line is the corresponding 
7—ray flux. The UHECR flux observed by HiRes 1] and two 
estimates of the extragalactic diffuse 7—ray background de- 
duced from EGRET data as blue (higher) [g] and red (lower) 
crosses Hi are also shown. 



The Diffuse GeV-TeV 7— ray flux. We now discuss the 
possible range of contributions of UHECR to the diffuse 
7—ray flux in the EGRET band. In Fig. Q] we show a 
scenario where the dip is due to pair production by ex- 
tragalactic protons. The lower red (solid) line was fitted 
to the HiRes spectrum [l[ at energies E > 2 EeV and the 
corresponding EM cascade flux is shown as blue (dash- 
dotted) line. Fig. Q] shows that practically all EM energy 
ends up in the GeV-TeV region. By also showing the pro- 
ton spectrum (green, dashed line) and the corresponding 
cascade flux (magenta, dotted line) when pair produc- 
tion by protons is neglected, Fig. Q] demonstrates that 
the GeV-TeV 7—ray flux is dominated by pair produc- 
tion losses of protons in this scenario. 

The flux in energy carried by a differential spectrum 



j(E) is given by / E 2 j(E)dlnE. One can then see from 
Fig. Q] that the UHECR energy lost to pair production 
by protons (energy flux difference between the red, solid 
and the green, dashed line) appears as the energy flux dif- 
ference between the blue, dash-dotted and the magenta, 
dotted line, up to the pair flux which is not shown. Fur- 
thermore, the UHECR energy lost to pion production by 
protons (energy flux difference between the two red, solid 
lines) appears as the energy flux in the magenta, dotted 
line, up to the neutrino flux which is not shown. 
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FIG. 2: Same as Fig. [T] but for a scenario with power law 
injection oc E~ 2 and evolution oc (1 + z) 3 . 

In Fig. [5] we show a scenario where the UHECR flux 
is extragalactic only above the ankle. In the particular 
case a = 2 one can see the relation between the energy 
deposited in GeV-TeV photons and the energy lost by 
protons more directly. For this purpose we artificially 
continue the proton flux to low energies where it has the 
same order of magnitude as the secondary photon flux. 
Since protons loose similar amounts of energy to pion 
and to pair production in this scenario, the contribution 
of these processes to the GeV-TeV 7—ray flux is also 
comparable, contrary to the dip scenario of Fig. [TJ 

In Fig. [3] we show the range of possible contributions of 
UHECR interactions to the EGRET flux, as a function 
of redshift evolution index to. We express this as the 
fraction of integral fluxes between 1 and 2 GeV, where 
the EGRET energy flux E 2 j(E) is minimal. Fig.[3]shows 
that the diffuse 7—ray flux in the EGRET band strongly 
depends on the source luminosity evolution index to. In 
the ankle scenario, for a given value of to, it still depends 
on other parameters within a factor 3, whereas in the 
dip scenario the scatter is much smaller due to partial 
degeneracy between to and a. All realistic astrophysical 
source distributions have m > 0, which implies that the 
contribution of secondary photons from UHECR will be 
at least ~1%. For stronger evolution, to = 3 — 4, this 
fraction increases to more than 50 % of the EGRET flux. 
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FIG. 3: Dependence on redshift evolution index m of the 
minimal and maximal fractional contribution of UHECR in- 
teractions to the EGRET flux between 1 and 2 GeV. Blue, 
dashed lines are for fitting above 2 EeV (dip scenario of the 
type shown in Fig. [lj and red, solid lines are for fitting the 
UHECR spectrum above 40 EeV. 
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FIG. 4: As Fig. but as function of the UHECR injection 
spectral index a. 



In Fig. 5] we show the range of possible contributions of 
UHECR interactions to the EGRET flux, as a function 
of the UHECR injection power law index a. Contrary 
to the case of Fig. [3l the scatter is larger, especially for 
small values of a = 2 — 2.4. This is due to the strong 
dependence of the flux in the EGRET band on the value 
of m for any given a, see Fig. |3] The lower lines corre- 
spond to minimal values of m = —2, while the maximum 
is defined by m — 4 for a < 2.4, and by smaller values 
of m for a > 2.4. Other parameter combinations would 
overproduce the cosmic ray flux below ~ 10 EeV. 

In Fig. [5] we compare the range of EM cascade fluxes 
from UHECR with other possible astrophysical contri- 
butions in the EGRET band. Note that most of the 
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FIG. 5: The possible range of UHECR induced cascade 
fluxes (light shaded band) compared to estimated 7— ray 
fluxes directly produced by starforming galaxies [6|, star- 
bursts 0] , large scale structure formation shocks || , AGNs [j| , 
and 7— ray bursts [ill ]. The dark shaded band shows the 
range of EM cascade fluxes from UHECR sources evolving as 
AGN [2l]]. The extragalactic diffuse 7— ray background from 
EGRET is as in Fig. Q] 

uncertainty of the UHECR cascade flux comes from the 
unknown source evolution. The scatter for a given evolu- 
tion such as for AGN is thus much smaller, as seen from 
Fig. [5j We remark that dark matter annihilations may 
also contribute to the diffuse flux [22I ]. 

The GeV-TeV cascade flux in scenarios where extra- 
galactic cosmic rays dominate down to below the an- 
kle at ~ 5 x 10 18 eV would be practically diffuse for an 
instrument such as GLAST which will be sensitive to 
anisotropies down to ~ 0.1% [2Cj. This is because the 
cosmic ray flux below the ankle is dominated by cosmo- 
logical sources and is isotropic at the percent level, and 
because large scale magnetic fields should lead to signifi- 
cant additional isotropization of primary protons [23| and 
secondary pairs. 

In contrast, lacking statistics, UHECR anisotropies at 
the 10% level can currently not be ruled out at ener- 
gies around or above the GZK cutoff, E > 40 EeV, be- 
cause sources of such UHECR may have a small den- 
sity ~ 10" 5 Mpc" 3 Q. Furthermore, at energies E ~ 
200 GeV-10 TeV the 7— ray absorption length in the IR 
background becomes small compared to the Hubble ra- 
dius, but still large compared to the ~ lOMpc length 
scales of pion production and EM cascade development 
around the source. This is also reflected in the amount 
of 7— ray flux suppression in Figs. [T] and [2j Independent 
of the poorly known size of deflection, this could lead 
to considerable correlation with nearby UHECR sources, 
as in the case where discrete sources emit very high en- 
ergy 7— rays directly [25|. In scenarios such as in Fig. [21 
where the secondary cascade flux is dominated by pion 
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production due to relatively hard UHECR injection spec- 
tra, this flux could, therefore, exhibit detectable small 
scale anisotropy around ~ 200 GeV. 

Conclusions. UHECR interactions with low energy 
photons can significantly contribute to the observed dif- 
fuse flux of 7— rays at energies between ~100 MeV and 
~ TeV. In this paper we studied the dependence of this 
contribution on unknown parameters of astrophysical 
UHECR scenarios. We found that UHECR contribute 
no less than 1% to the observed EGRET flux, and up 
to 50% in some cases. This suggests that the GLAST 
satellite, which at GeV energies will be ~ 30 times more 
sensitive to point sources than the EGRET experiment, 
will likely be sensitive to the UHECR induced contribu- 
tion. Even ground-based instruments such as HESS and 
the future CTA may be sensitive to the cascade flux be- 
tween ~ 0.1 and ~ 10 TeV, although such experiments 
are less well suited for diffuse backgrounds. 

If sources of extragalactic highest energy cosmic rays 
are rare and dominate the flux down to only ~ 5 x 
10 18 eV, the cascade background may have significant 
anisotropy at energies around 200 GeV. 

To summarize, future measurements of resolved and 
unresolved components of the diffuse EGRET 7— ray 
background or upper limits on such components can give 
important information on UHECR origin and the distri- 
bution of their sources. 
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